clc;
clear;

pde.domain = [0 1 0 1 0 1]
h = 1
bt = ["P1" "P2"]

for i = bt
% mesh = process_mesh(pde, h);

% fs = process_fespace(mesh, i);

% process_showmesh(mesh, fs, [1 1 1 1]);

process_showbasis(i, "x")
end
